F/6  3/1 


AD-A074  640 


AEROSPACE  CORP  EL  SE6UND0  CA 

EQUATORIAL  INFRARED  CATAL06UE #  NUMBER  1»(U) 

JUN  78  L  H  SWEENEY t  T  F  HElNSHEIMER 


F04701-C -77-0078 


UNCLASSIFIED  TR-0078C 3409-20 )-l  6IDEP-E135-0405 

i  of  i  B^egpvflB 

, . — flj  ~aMI  Eeas^easi!  P  -'-Ui  Ii  jjjl  I  - 5  1 


NL 


DOC  FILE  COPY  #JA0  74640 


/ 


7  T 
J  i 


>£Pa 


tz 


/ 


INTERIM 


Aerospace  Report  No. 
,  UJ  TR-00  78(3409 -2jD)-l  ( 
Formerly 

TOR -00  78(3409- 20)- 2 


EQUATORIAL  INFRARED  CATALOGUE,  NUMBER  1 


Prepared  by 


L. 


H.  Sweeney,  T.  F.  Heinsheimer,  and  F.  F.  Yates^, 
The  Aerospace  Corporation 
and 

<S>  P./Maran  R.,  Lesh/  and  T.  A. /Nagy 

/  yU. . 


j/  15  Jun«Tiit978 


Programs  Group 

THE  AEROSPACE  CORPORATION 
El  Segundo,  Calif.  90245 


Prepared  for 

SPACE  AND  MISSILE  SYSTEMS  ORGANIZATION 
AIR  FORCE  SYSTEMS  COMMAND 
LOS  ANGELES  AIR  FORCE  STATION 
P.  O.  Box  92960,  Worldway  Postal  Center 
Los  Angeles,  California  90009 


Contract  No.  'F04701  -C-77-0078 

J  : 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


D  D  C 

OCT  4  1979 

nstolhinrtsili 

A 


1 


GOVERNMENT. INDUSTRY  DATA  EXCHANGE  PROGRAM 


GENERAL  DOCUMENT  SUMMARY  SHEET 


1  OF  2_ 


I.  ACCESS  NUMBER 


EI3S-0<40S 


Pitas*  Typo  All  Information  —  Soo  Inatnjdions  on  Rovorso 


2.  COMPONENT/ PART  NAME  PER  GIDEP  SUBJECT  THESAURUS 

Environmental  Studies,  Screening,  Statistical 


3.  APPLICATION 


4.  MFR  NOTIFICATION 


IS.  DOCUMENT  ISSUE  month/ Yoar) 


Engineering 


6.  ORIGINATOR'S  DOCUMENT  TITLE 


|"xl  not  applicable  June  1978 


Interim  Equatorial  Infrared  Catalogue,  Number  1 


8.  ORIGINATOR  S  DOCUMENT  NUMBER 

TR-0078 ( 3^09-20 )-l 


10.  DOCUMENT  (SUPERSEDES)  (SUPPLEMENTS)  ACCESS  NO. 


9.  ORIGINATOR'S  PART  NAME/IDENTIFICATION 


12.  MANUFACTURER 


113.  MANUFACTURER  PART  NUMBER 


7.  DOCUMENT  TYPE 


NONSTD  PART  SPEC 


114.  industry/government  standard  number 


r  CK.  i'f,  I©  A  laJC 


^5.  OUTLINE,  TABLE  OF  CONTENTS,  SUMMARY,  OR  EQUIVALENT  DESCRIPTION 

>  \ 

\  S  A.  ^  I  O  ^  c-vJ £  \ 

''“^Provides  a  list  of  Infrared  sources,  detected  at  wavelength  2.7  micrometer. 
Confidence  is  achieved  by  multiple  observations  of  each  listed  source  and 
by  careful  computer  processing  of  the  satellite  observational  data.  Every 
source  has  been  observed  on  two  or  more  days  or  at  times  on  the  same  day, 
at  least  20  min  apart.  Generally,  observed  positions  are  precessed  to  a 
common  year  and  averaged  whenever  two  or  more  measurements  are  within  U5 
arc  sec  of  one  another.  The  resulting  number  of  apparent  infrared  positions 
exceeds  5000.  The  halfwidths  of  U5  arc  sec  for  positional  averaging  is  rath¬ 
er  large  in  order  to  simplify  the  data  processsing  logic. __ 

'\ 


16.  KEY  WORDS  FOR  INDEXING 


Infrared  Sources;  Astronomy 


(Doc  Des— P) 


17.  GIDEP  REPRESENTATIVE 


■  18.  PARTICIPANT  ACTIVITY  AND  CODE 


James  M.  Teresi 


Aerospace  Corp.,  Los  Angeles,  CA  (BA) 


FORM 
1  NOV  7* 


REPRODUCTION  OR  DISPLAY  OF  THIS  MATERIAL  FOR 
SALES  OR  PUBLICITY  PURPOSES  IS  PROHIBITED 


INSTRUCTIONS  FOR  COMPLETING  THE  GIDEP  GENERAL  DOCUMENT  SUMMARY  SHEET 

NOTE:  Completion  of  a  Summary  Sheet  by  the  participant  is  not  mandatory  for  document  acceptance  into  GIDEP.  A  Summary  Sheet  will  be  prepared  by  the 
GIDEP  Operations  Center  for  document  submittals  received. 

PAGE  1  OF  Enter  the  total  number  of  summary  sheet  pages 


BOX 

1  Leave  blank — entry  will  be  completed  by  GIDEP  Operations  Center. 

2  Enter  standard  nomenclature  associated  with  GIDEP  Subject  Thesaurus  selected  from  Section  12,  Policies  and  Procedures  Manual. 

3  Indicate  application  which  the  device  was  used  (e.g.,  ground,  missile,  shipboard,  spacecraft;  refer  to  P  &  P  Manual,  Section  13). 

4  Device  manufacturer  must  be  notified  of  test  results.  Manufacturer  approval  of  report  is  not  required — include  pertinent  manufacturer 

correspondence  with  document  submittal  to  GIDEP;  check  NOTIFIED  entry.  If  document  is  for  a  nonstandard  part  or  of  a  general 
nature  and  a  manufacturer  is  not  identified,  check  NOT  APPLICABLE. 


5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


16 


17 

18 


Enter  month  and  year  of  document  issue. 

Enter  complete  document  title  exactly  as  it  appears  on  originator  document. 

Identify  document  type  by  inserting  letter  X  by  appropriate  descriptor. 

Enter  document  number  exactly  os  it  appears  on  originotor  document. 

Enter  part  name  and  identification  as  assigned  by  organization/agency  originating  the  report;  if  not  specified,  enter  N/A  (Not  Applicable). 

Delete  either  SUPERSEDES  or  SUPPLEMENTS.  If  document  supersedes/supplements  on  existing  GIDEP  document,  enter  GIDEP  microfilm 
access  number  of  appropriate  document.  If  document  neither  supersedes  nor  supplements  an  existing  GIDEP  document,  enter  the  word  NONE. 

Enter  the  single  symbol  coding  for  environmental  exposure  as  defined  in  Subject  Thesaurus,  Section  12,  P&P  Manual  (e.g.,  C  -  Salt  Spray; 

V  -  Vibration;  %  -  Shelf  Life);  if  not  specified,  enter  N/A  (Not  Applicable). 

Enter  manufacturer  abbreviation  and  H-4  Code  number  listed  in  GIDEP  Manufacturer  List.  If  manufacturer  is  not  listed,  enter  the  phrase, 

SEE  BOX  15,  enter  manufacturer’s  full  name  and  division  (if  any)  in  Box  15.  If  more  than  one  manufacturer,  enter  phrase,  SEE  BOX  15;  enter 
additional  manufacturers  as  appropriate.  If  manufacturer  is  not  specified,  enter  N/A  (Not  Applicable). 

Enter  complete  part  number.  Use  open  0  for  alpha  letter  0,  and  use  0  for  numeric  zero.  If  more  than  one  part  number,  enter  phrase,  SEE  BOX 
15;  enter  additional  part  number(s).  If  a  part  number  is  not  specified,  enter  four  dashes  ( - ). 

Enter  standard  part  number  such  as  the  IN  or  2N — diode  and  transistor  designators.  For  GIDEP  purposes,  any  military  assigned  number  is 
considered  as  a  government  standard  part  number.  If  more  than  one  standard  number,  enter  phrase  SEE  BOX  15;  enter  additional  standard  number(s). 

If  subject  matter  in  document  can  be  categorized  into  more  than  one  major  subject  category,  enter  additional  subject  categories  in  upper  right- 
hand  corner.  Briefly  summarize  test  results  or  material  detailed  in  text  of  document.  Include  any  pertinent  details  or  comments  required  for 
proper  interpretation  of  material  presented  (e.g.,  peculiar  environmental  capabilities,  unique  electrical  characteristics  that  may  be  "state-of-the- 
art,"  or  characteristics  that  restrict  part  usage  to  particular  applications  or  any  other  details  that  may  aid  a  prospective  user  of  the  part). 

Enter  appropriate  words  or  phrases  that  enhance  information  retrieval  on  subject  matter(s)  contained  in  document.  As  a  secondary  data  retrieval 
technique  within  each  applicable  Major  Category  (entry  2),  the  document  is  referenced  in  the  computer  doto  bonk  ond  Report  Index  according  to 
each  key  word.  Do  not  use  abbreviations  or  words  that  are  part  of  the  subject  category  listed  in  Box  2.  Key  word  phrases  are  limited  to  60  total 
characters  and  blank  spaces.  Separate  key  words  and/or  phrases  with  commas. 

Enter  signature  or  name  of  GIDEP  Representative. 

Enter  name,  city,  and  state  of  participant  activity  or  corporation  and  division  submitting  the  document  and  GIDEP  two-character  code  (e.g.,  XI). 


OD  FORM  2000  (t  NOV  7«)  (RACK) 


'T~ 


CONTENTS 


FOREWORD .  iii 

1.  INTRODUCTION .  1 

2.  GENERAL  DESCRIPTION  .  2 

2.  I  EIC  Processing  Steps .  9 

2.1.1  Step  1  of  3  . , .  11 

2.  1.2  Step  2  of  3  .  13 

2.1.3  Step  3  of  3 .  17 

3.  SUMMARY .  19 

REFERENCES  .  20 

COLUMN  HEADINGS  FOR  EIC-1 .  21 


iv 


TABLES 


1.  EIC-1  Parameters .  3 

2.  Satellite  Observations  by  Date  and  Declination  Zone .  4 

3.  EIC-1  Sources  by  Declination  Zone  and  Flux 

Density  Interval  .  6 

4a.  Crosscorrelations  for  137  EIC-1  Stars  at 

Declination  +8  deg . 7 

4b.  EIC-1  Crossreference  and  Coordinate  Difference 

Statistics .  8 

5.  Identified  EIC-1  Sources  by  Spectral  Type  and  Flux 

Density  Interval .  10 

6.  EIC-1  Flux  Density  by  Galactic  Latitude  Zone .  10 

7.  EIC-1 .  24 

FIGURES 

1.  EIC  Processing  -  Step  1  of  3 .  12 

2.  EIC  Processing  -  Step  2  of  3 .  14 

3.  Typical  Alignment  Angular  Variation  .  16 

4.  EIC  Processing  -  Step  3  of  3 .  18 


v 


1.  INTRODUCTION 


The  Equatorial  Infrared  Catalogue,  Number  1  (EIC-1,  pronounced 

ICE-ONE)  lists  positions  for  896  sources  detected  at  wavelength  2.7-pm 

with  U.S.  Air  Force  satellite  sensors  within  10  deg  of  the  celestial  equator. 

EIC-1  together  with  planned  updates  will  provide  a  systematic  survey  of 
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sources  with  flux  densities  greater  than  4  X  10  W  cm  jxm  in  this 
equatorial  zone.  The  satellite  data  are  processed  to  provide  flux  density  to 
a  root  mean  square  (rms)  precision  of  14  percent  and  positional  coordinates 
to  a  rms  accuracy  of  2  arc  sec.  Positions  listed  in  EIC-1  are  computed 
relative  to  those  of  the  Astronomische  Gesellschaft  Katalog  3  (AGK3)  (Ref.  1) 
and  the  Smithsonian  Astrophysical  Observatory  (SAO)  Star  Catalogue  (Ref.  2). 
These  catalogues  use  equinox  1950.0  and  the  coordinate  system  of  The  Fourth 
Fundamental  (FK4)  Catalogue  (Ref.  3). 

The  systematic  positional  errors  of  EIC-1  with  respect  to  a  sample 
of  658  reference  stars  (500  from  the  AGK3  and  158  from  the  SAO)  are  less 


2.  GENERAL  DESCRIPTION 


EIC-1,  the  first  in  a  series  of  2.  7-(xm  catalogues,  provides  a  list  of 
infrared  sources  that  is  assembled  with  high  confidence.  Considerable  effort 
has  been  expended  to  deal  with  the  known  problem  of  spurious  sources  in  space 
infrared  surveys.  Confidence  is  achieved  by  multiple  observations  of  each 
listed  source  and  by  careful  computer  processing  of  the  satellite  observa¬ 
tional  data.  Every  EIC-1  source  has  been  observed  on  two  or  more  days  or 
on  the  same  day  at  times  at  least  20  min  apart.  NASA-GSFC  correlation  of 
EIC  sources  with  photographic  plate  images  and  multichannel  photometry  of 
selected  sources  with  a  ground-based  infrared  telescope  will  be  addressed  in 
a  later  publication.  Selected  results  of  the  NASA  work  appear  in  References 
4  and  5.  EIC-1  parameters  are  given  in  Table  1. 

The  development  of  EIC-1  began  with  the  processing  of  satellite  data 
recorded  in  November  1976.  The  sky  coverage  at  that  time  emphasized  the 
+8-deg  declination  band.  During  the  14.  2  hr  recording  interval,  1.  3  hr  on  the 
9th  and  12.9  on  the  10th,  more  than  2400  of  some  3400  measurements  were 
confined  to  one  deg  of  declination  between  +8  and  +9  deg. 

Subsequent  recordings  still  emphasize  northern  declinations,  which 
means  that  the  survey  will  be  complete  in  these  zones  prior  to  the  others. 

Table  2  shows  by  Greenwich  date  and  declination  zone  the  number  of  measure¬ 
ments  processed  to  make  EIC-1.  Generally,  observed  positions  are  precessed 
to  a  common  year  and  averaged  whenever  two  or  more  measurements  are 
within  45  arc  sec  of  one  another.  The  resulting  number  of  apparent  infrared 
positions  exceeds  5000.  The  halfwidths  of  45  arc  sec  for  positional  averaging 
is  rather  large  in  order  to  simplify  the  data  processing  logic.  Large  positional 
variances  are  investigated  to  see  that  close  sources  are  not  merged,  but  so 
far  no  sources  within  100  arc  sec  of  each  other  have  been  seen.  Most  of  the 
apparent  source  positions  have  not  been  redetected  and  hence  are  omitted 
from  EIC-1.  Omitted  sources  either  are  spurious  or  require  confirmation  by 
further  observation. 


Table  2.  Satellite  Observations  by  Date  and  Declination  Zone 


The  lack  of  uniformity  in  sky  coverage  is  easily  seen  in  Table  3, 
which  shows  the  count  of  sources  by  3-deg  wide  declination  zones  in  two  flux 
density  intervals.  For  comparison,  the  counts  of  2.  2-pm  sources  are  shown 
from  the  Two-Micron  Sky  Survey  (generally  abbreviated  as  TMSS)  (Ref.  6). 

If  the  sky  coverage  were  uniform,  approximately  equal  numbers  of 
EIC-1  sources  would  be  found  in  the  respective  declination  bands.  Every  real 
TMSS  source  within  the  declination  range  of  the  EIC  survey  should  eventually 
be  included  in  the  EIC.  Table  3  provides  data  on  the  percentage  of  TMSS 
sources  (by  declination  band)  that  have  been  observed  with  the  satellite  sensors 
and  that  meet  the  criteria  for  inclusion  in  EIC-1.  The  deviations  of  these  num¬ 
bers  from  100  percent  (neglecting  the  small  number  of  possibly  spurious 
sources  in  the  TMSS)  constitute  an  indication  of  the  lack  of  uniformity  in  sky 
coverage  of  EIC-1.  An  EIC-1  identification  with  a  TMSS  source  requires 
positional  coincidence  within  180  arc  sec;  where  such  identifications  are  valid 
the  EIC-1  position  constitutes  a  preferred  estimate  that  can  be  used  for  further 
ground-based  observations  and  optical  and  radio  source  identification  purposes. 
The  number  of  EIC-1  sources  with  positions  that  match  SAO  catalogue  and 
AGK3  positions  to  within  18  arc  sec,  the  number  that  match  AFGL  (Ref.  7) 
positions  to  within  360  arc  sec,  and  the  number  that  match  TMSS  positions 
with  180  arc  sec,  are  summarized  in  Tables  4a  and  4b. 

Table  4a  concentrates  on  the  +8-deg  declination  zone  where  EIC-1  is 
only  76  percent  complete  to  the  TMSS  flux  threshold,  since  just  60  of  the  79 
TMSS  stars  were  matched.  These  60  represent  44  percent  of  the  EIC-1  stars 
in  the  +8-deg  zone.  The  ragged  sky  coverage  will  become  more  and  more 
uniform  as  additional  satellite  data  are  received  and  processed.  Table  4b 
shows  crosscorrelations  and  coordinate  difference  statistics  for  all  EIC-1 
star  matches  within  10  deg  of  the  celestial  equator.  The  number  of  EIC-1 
sources  that  do  not  match  those  of  any  star  in  these  four  catalogues  is  88. 

These  deserve  fuither  study;  preliminary  observations  by  S.  P.  Maran,  S.  G. 
Kleinmann  and  R.  Joyce  show  that  some  of  them  have  interesting  infrared 
spectra. 
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Table  3.  EIC  - 1  Sources  by  Declination  Zone  and 
Flux  Density  Internal 


Flux  Density 
10-16\y  cm_2  pim "  1 

Declination 
Epoch  1950.0 

3  to  74 

75  and  More 

EIC- 1 
Count 

TMSS 

Count 

EIC  - 1 
Coverage 
of  TMSS  (%) 

Above  +9° 

28 

3 

31 

* 

* 

+6°  to  +9° 

309 

31 

340 

175 

87 

+3°  to  +6° 

186 

23 

209 

164 

79 

0°  to  +3° 

67 

15 

82 

155 

40 

» 

LO 

o 

r-t- 

o 

o 

o 

33 

37 

70 

161 

37 

-6°  to  -3° 

42 

17 

59 

166 

29 

-9°  to  -6° 

66 

33 

99 

175 

43 

Below  -9° 

4 

2 

6 

>;c 

* 

T  o tal s 

735 

161 

896 

Not  of  interest  in  present  context 


Table  4a.  Crosscorrelations  for  137  EIC-1  Stars 
at  Declination  +8  deg 


Reference 

Catalogue 

(Wavelength) 

Reference 

Stars 

Matches 

Match 
Criteria 
(arc  sec ) 

AGK3  (visual) 

3284 

89 

18 

SAO  (visual) 

2073 

83 

18 

TMSS  (2.2  |rm) 

79 

60 

180 

AFGLi  (4.  2  nun) 

26 

20 

360 

19.  4 

rsj 

sO 

r- 

ro 

00 

o 

• 

• 

• 

• 

ro 

fM 

00 

03 

ro 

00 

• 

• 

• 

O 

i 

o 

i 

o 

1 

CO 

03 

■«■ 

•3E- 

1 

u 

X 

o 

lO 

in 

(J 

o 

CO 

w 

ctJ 

S 

LO 

NO 

m 

(M 

<v 

u 

C 

CO 

d> 

D 

4-> 

'-*-4 

V 

cc 

o 

c 

• 

o 

in 

•  H 

4J 

rd 

v  o 

c 

d  ^ 

•  H 

o  ^ 

i—4 

u 

N  -g 

4) 

o 

A 

a 

U 

03 

o 

Tt< 

O 

r~ 

OJ 

o 

oo 

o 

•t 

•» 

oo 

00 

ro 

rO 

o 

o 

O 

o 

O 

O 

O 

o 

■*H 

■*— t 

4~ 

-f 

4- 

4- 

o 

0 

O 

O 

4-> 

4-» 

4-» 

o 

o 

O 

O 

03 

o 

O 

o 

1 

^rH 

1 

^■H 

1 

■H 

1 

*-i  m  r 

M  .rH  1> 
"  h  TJ 

II  II  II 


3-1  stars 

AFGL 

1010 

1809 

1825 

HH 

W 

o 

w 

cn 

O' 

i 

£ 

s 

■»-«  1  I  I 

U-, 

H 

o 

0 

+ 

CO 

4> 

00  m  (M  0s 
M  in  vo  r- 

w  "-1  ^  'l'  't 

^  o&  o&  o<3  o<5 

t-J  r-  ro  oo  oo 

<M  in  in  r- 

■H 


633  &  635  +10353 

659  &  660  -10424  2164 

678  &  679  -  2195 


Once  a  position  is  matched,  additional  data  from  the  appropriate 
reference  catalogue  matched  are  appended  to  the  EIC  data  record.  Of  par¬ 
ticular  interest  in  the  associated  information  are  the  stellar  spectral  type 
and  proper  motion.  The  spectral  type  is  obtained  from  the  Henry  Draper 
(HD)  Catalogue  (Ref.  8)  using  the  HD-SAO-DM  Cross  Index  (Ref.  9)  whenever 
an  EIC  star  is  identified  with  an  SAO  star  in  the  Cross  Index.  If  not  found 
there,  the  spectral  type  is  taken  from  the  AGK3,  or  as  a  last  choice,  is  from 
the  SAO  without  use  of  the  Cross  Index.  Table  5  displays  the  count  of  sources 
by  spectral  type  in  two  flux  density  intervals.  As  expected  for  the  2.  7-pm 
EIC  wavelength,  nearly  all  the  sources  are  associated  with  type  K  and  M  stars. 

Proper  motion  is  obtained  from  the  AGK3  as  first  choice  and  the  SAO 
as  last  choice.  The  proper  motion  is  not  significant  for  most  of  the  star  ob¬ 
served.  However,  it  cannot  be  ignored  completely.  The  current  technique  of 
precessing  all  positions  to  a  common  equinox  and  then  averaging  the  coordinates 
may  be  modified  in  the  future  to  account  for  proper  motion  as  the  span  of  time 
for  EIC  recordings  grows,  since  several  of  the  stars  observed  do  move  more 
than  1  arc  sec  per  year.  The  proper  motion  for  stars  is  always  used  to  deter¬ 
mine  infrared  sensor  alignment  (discussed  below). 

Galactic  coordinates  are  computed  and  listed  for  each  source  position. 
Table  6  shows  EIC-1  flux  density  distribution  by  galactic  latitude  zone.  As 
expected  for  real  sources  in  the  Galaxy,  their  number  decreases  as  the  dis¬ 
tance  from  the  galactic  equator  increases. 

2.  1  EIC  PROCESSING  STEPS 

The  EIC-1  listing  is  produced  in  three  distinct  processing  steps.  The 
first  step  is  repeated  for  each  digital  tape  of  satellite  measurements.  After 
all  the  tapes  in  a  continuous  series  are  processed  through  step  1  or  the  end 
of  an  orbital  cycle  of  observations  is  reached,  whichever  occurs  first,  step  2 
is  performed  for  the  interval  just  defined.  The  accumulated  outputs  from  all 
of  the  step  2  processing  constitute  the  input  to  step  3.  In  step  3,  statistics  are 
computed,  sources  are  screened,  and  catalogue  numbers  are  assigned. 
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Table  5.  Identified  EIC-1  Sources  by  Spectral  Type 
and  Flux  Density  Interval 


Flux  Density 

. n- 16  -2  -1 
10  W  cm  pm 

Spectral  Type 

Type 

Not 

Listed 

T  otals 

O 

B 

A 

F 

G 

K 

M 

R 

N 

S 

75  and  above 

0 

1 

0 

1 

2 

28 

97 

0 

5 

1 

26 

161 

3  to  74 

0 

1 

6 

8 

26 

300 

177 

1 

2 

0 

214 

- 

735 

T  otals 

0 

2 

6 

9 

28 

328 

274 

1 

7 

1 

240 

896 

i 

i 


Table  6.  EIC-I  Flux  Density  by  Galactic  Latitude  Zone 


2.1.1  Step  1  of  3 

The  survey  is  being  carried  out  using  infrared  and  visual  wavelength 
sensors  aboard  U.  S.  Air  Force  satellites.  Data  processing  of  a  satellite's 
visual  wavelength  sensor  determines  the  orientation  of  that  sensor  as  a  func¬ 
tion  of  time  where  the  true  positions  of  the  visual  stars  observed  are  obtained 
from  the  SAO  star  catalogue  (see  Figure  1).  Approximately  1000  of  the 
brighter  SAO  stars  distributed  throughout  the  sky  provide  the  reference  posi¬ 
tions  to  estimate  the  visual  sensor  orientation.  Digital  tapes  supply  infrared 
sensor  data,  visual  sensor  orientation  angles,  and  satellite  location.  These 
input  tapes  are  processed  chronologically  until  the  end  of  an  observation 
interval  or  a  satellite  orbital  cycle  is  completed.  This  processing,  referred 
to  as  step  1,  removes  sensor  design  features,  detects  apparent  celestial 
sources,  and  determines  the  1975.0  position  for  each  detected  source. 

The  satellite's  infrared  sensor  is  a  linear  array  of  detectors  stag¬ 
gered  in  position  to  provide  overlapping  coverage  as  the  sensor  scans  across 
the  sky.  This  configuration  produces  multiple  detections  that  must  be  com¬ 
bined  to  form  representative  measurements  as  adjacent  detectors  may  re¬ 
spond  to  only  a  fraction  of  a  source  image.  A  representative  measurement  is 
computed  from  up  to  three  single  detector  responses  that  are  closely  spaced 
in  time  and  position.  Representative  time  is  the  time  of  the  maximum  flux 
measurement,  representative  flux  density  is  the  maximum  flux  measurement, 
and  representative  position  is  computed  as  a  weighted  average  using  flux 
density  weighting  and  the  midpoints  of  responding  detectors.  Of  the  45,617 
representative  measurements  available  for  EIC-1,  only  29,285  are  associated 
with  sources  listed  in  EIC-1. 

An  apparent  celestial  source  is  detected  as  follows.  Detector  re¬ 
sponses  collected  over  a  span  of  several  minutes  of  sky  observation  are 
processed  as  a  group  where  each  response  is  checked  for  neighbors  in 
position  and  time.  An  apparent  celestial  source  detection  is  achieved  when 
three  responses  are  found 'within  a  40  by  40  arc  sec  sky  area  and  the  time 
between  the  first  and  last  of  these  responses  is  no  more  than  32  sec.  All 
responses  in  the  group  being  processed  that  fall  within  the  sky  area  of  the 
detection  are  then  combined  as  appropriate  to  form  representative  measure - 
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Figure  i.  EIC  Processing  -  Step  1  of  3 


ments.  As  many  as  three  responses  can  be  combined  to  form  a  single  re¬ 
presentative  measurement,  however,  most  measurements  are  computed 
from  just  one  or  two  responses.  The  29,285  representative  measurements 
associated  with  EIC-1  sources  resulted  from  47,868  detector  responses: 

12,404  singles,  13,  595  doubles,  and  only  2758  triples. 

Observed  positions  are  transformed  from  sensor  coordinates  to 
celestial  coordinates  using  an  approximate  transformation  between  the  visual 
and  infrared  reference  systems.  This  approximation  is  removed  in  step  2. 

Correction  for  annual  aberration,  nutation,  and  precession  from  the 
day  of  observation  to  the  nearest  year  (e.  g. ,  to  1977.0)  is  made  using  Inde¬ 
pendent  Day  Numbers  as  listed  in  the  American  Ephemeris  and  Nautical 
Almanac  for  the  appropriate  year  and  day.  A  correction  for  precession  from 
the  nearest  year  to  1975.0  coordinates  is  then  made  using  the  equations  and 
constants  of  The  Explanatory  Supplement  to  the  Astronomical  Ephemeris  and 
American  Ephemeris  and  Nautical  Almanac.  These  same  equations  and  con¬ 
stants  are  used  to  precess  reference  star  positions  from  1950.0  to  1975.0  for 
sensor  alignment  estimation  in  step  2.  Representative  measurements  pro¬ 
duced  by  successive  executions  of  step  1  processing  are  merged,  then  ordered 
by  right  ascension  in  preparation  for  step  2. 

2.1.2  Step  2  of  3 

Step  2  processing  (see  Figure  2)  corrects  for  two  different  kinds  of 
positional  errors.  The  first  of  these  errors  and  smaller  of  the  two  is  the 
correction  for  aberration  due  to  satellite  orbital  velocity  ("orbital  aberration"). 
It  is  a  function  of  the  satellite's  position  and  velocity  with  respect  to  the  center 
of  the  earth.  More  important  is  the  correction  for  sensor  alignment,  which  is 
required  due  to  variations  in  the  satellite's  temperature.  Temperature  vari¬ 
ations  cause  a  slow  change  in  the  optical  alignment  of  one  sensor  with  respect 
to  another.  The  alignment  drift  is  estimated  by  determining  values  for  two 
angles  defining  the  infrared  sensor  positional  deviation  from  a  satellite  refer¬ 
ence  frame  whose  orientation  in  space  is  known  by  way  of  the  visual  star 
sightings. 


A  typical  plot  of  infrared  sensor  angular  deviation  from  its  satellite 
reference  frame  over  a  complete  orbital  cycle  is  given  in  Figure  3,  Correc¬ 
tions  to  these  angles  are  estimated  for  each  orbit  since  changes  greater  than 
1  arc  sec  can  occur.  These  angles  and  the  corrections  are  very  small  so  that 
no  significant  loss  of  precision  is  noticed  when  small  angle  approximations 
are  used  in  the  transformation  of  errors  in  right  ascension  and  declination  to 
those  in  the  satellite  reference  frame,  or  the  use  of  each  angle  in  radians  to 
represent  its  sine  and  a  value  of  unity  for  its  cosine.  These  last  approxima¬ 
tions  are  used  in  the  equation  seen  in  Figure  3. 

The  selection  of  reference  stars  for  this  processing  was  done  iter¬ 
atively.  Any  EIC  star  with  multiple  observations  is  used  if  it  is  found  in 
either  the  AGK3  or  SAO  catalogue  and  it  is  not  too  bright  or  too  faint  (see 
below).  To  initialize  the  process,  SAO  stars  of  spectral  type  M  in  the  EIC 
declination  range  were  used.  Bright  stars  cause  extraneous  sensor  responses 
and  faint  stars  are  too  easily  associated  by  chance  with  thermally  induced  noise 
fluctuations  in  the  infrared  sensors.  The  growing  set  of  reference  stars 

currently  is  composed  of  615  stars  in  the  flux  range  of  6  to  6000  X  10 
-  2  -  1 

cm  pm  ,  470  AGK3  stars  and  145  SAO  stars  (434  of  the  AGK3  stars  are 
also  SAO  stars).  When  a  star  is  found  in  both  catalogues,  AGK3  positional 
coordinates  are  used.  EIC- 1  sources  with  fluxes  outside  the  above  range 
are  not  suitable  for  use  as  reference  stars. 

Observed  positional  errors  are  averaged  in  12-min  intervals  spanning 
an  orbital  cycle  and  then  smoothed  using  a  five-point  moving  average  before 
application  to  correct  the  observed  positions.  However,  any  interval  with 
fewer  than  five  measurements  is  treated  as  if  no  measurements  were  made. 
Intervals  with  no  measurements  are  supplied  values  using  linear  interpola¬ 
tion  between  surrounding  intervals  with  sufficient  measurements. 

It  should  be  noted  that  proper  motion  must  be  accounted  for  in  this 
processing  step  because  184  of  the  reference  stars  would  have  positional 
errors  greater  them  one  sec  of  arc  if  it  were  ignored  and,  consequently,  the 
sensor  alignment  correction  would  be  degraded.  Proper  motion  is  added  to  each 


15 


Typical  Alignment  Angular  Variation 


SAO  star  from  19  50.  0  to  the  date  of  the  satellite  observation  and  to  each  AGK3 
star  from  its  listed  date  of  observation  to  the  date  of  satellite  observation. 
Precession  to  epoch  1975.0  has  already  been  accounted  for  at  this  point  in  the 
processing.  Step  2  processing  outputs  are  accumulated  to  form  the  input  to 
step  3. 

2.  1.  3  Step  3  of  3 

Observations  input  to  step  3  are  ordered  to  facilitate  the  computation 
of  statistics.  The  observations  are  ordered  by  right  ascension  first,  then 
reordered  locally  by  declination  whenever  two  or  more  stars  have  observations 
interlaced  with  respect  to  right  ascension  alone.  The  ordered  observations 
are  then  processed  serially  to  determine  various  statistics:  average  position 
and  flux,  standard  deviations,  and  flux  summaries. 

Next,  high  confidence  sources  are  screened  (see  Figure  4).  Stars 
seen  on  two  or  more  days  or  twice  on  the  same  day  at  times  20  min  or  more 
apart  are  considered  for  further  processing.  A  measurement  not  associated 
with  a  multiply  observed  source  is  either  marked  as  bad  or  not  marked  at  all, 
but  held  for  subsequent  processing  together  with  future  observations.  Bad 
data  must  be  recognized  and  flagged.  Some  satellites  have  two  linear  arrays 
of  detectors  with  overlapping  sensitivity  ranges  and  when  a  star  is  seen  that 
exceeds  either  sensitivity  range,  bad  data  are  produced.  A  few  stars  exceed 
the  limits  of  both  arrays.  Bright  stars  produce  a  fringe  of  low-level  spurious 
detections  that  must  be  recognized  and  eliminated.  Each  flux  density  measure¬ 
ment  with  a  value  less  than  60  percent  of  the  average  flux  observed  for  a  source 
on  a  given  date,  or  less  than  20  percent  of  the  overall  average  flux,  is  assumed 
to  be  associated  with  this  low-level  fringe  or  other  extraneous  responses.  The 
average  flux  for  each  observational  interval  is  computed  sequentially  during 
this  processing.  This  logic  eliminated  nearly  10  percent  of  the  representative 
measurements.  The  moon  and  Venus  were  observed,  but  neither  of  them  is 
found  at  the  same  positions  in  the  sky  at  times  20  min  or  more  apart  in  the 
satellite  data,  nor  are  they  of  interest  in  the  present  context,  so  they  do  not 
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Figure  4.  EIC  Processing  -  Step  3  of  3 


appear  in  EIC-1.  Finally,  comparisons  are  made  to  star  positions  listed  in 
other  catalogues  and  the  EIC-1  star  numbers  are  assigned.  Star  positions 
are  compared  in  1975.0  coordinates  including  a  correction  for  proper  motion 
whenever  the  EIC-1  star  has  been  identified  with  an  AGK3  or  SAO  star. 

Proper  motion  is  not  used  in  the  calculation  of  1950.  0  positions  listed  for  the 
stars  in  EIC-1.  Catalogue  star  numbers  are  assigned  after  ordering  the  stars 
by  right  ascension  in  19  50.  0  coordinates. 

3.  SUMMARY 

The  completeness  of  this  survey  to  a  given  flux  density  level  is  not 
known.  Although  the  satellite  sensors  were  known  to  be  operating  w'hile  point¬ 
ing  at  the  sky  and  a  fundamental  threshold  is  known  for  each  sensor,  insuffi¬ 
cient  information  is  provided  in  the  recorded  data  received  for  EIC  process¬ 
ing  to  explain  why  certain  expected  measurements  are  not  present.  Complex 
satellite  and  ground  station  processing  algorithms  operate  on  the  sensor  data 
prior  to  it  being  recorded.  Sensor  data  is  lost  to  the  EIC  processing  when¬ 
ever  adaptive  threshold  logic  raises  thresholds  while  controlling  the  data  rate 
and  also  when  the  ground  processing  is  too  busy  to  record  sensor  pointing  data. 
About  once  every  30  min,  a  25 -sec  block  of  pointing  data  is  found  to  be  miss¬ 
ing  on  the  data  tapes. 

In  EIC-1  a  measure  of  the  completeness  to  the  threshold  level  of  the 
TMSS  was  made  by  comparing  positional  coordinates.  For  EIC -2  an  additional 
measure  will  be  added.  The  sky  within  10  deg  of  the  celestial  equator  will  be 
divided  into  1-deg  squares  and  the  minimum  recorded  flux  level  for  each 
square  will  be  displayed.  EIC -2  will  also  include  the  results  of  processing 
satellite  data  recorded  back  in  December  of  1975  as  well  as  some  recorded 
in  March  and  June  of  19  78,  all  in  addition  to  the  reprocessing  of  EIC-1  data. 
The  experience  gained  in  the  production  of  EIC-1  will  be  used  to  adjust  param¬ 
eter  values  in  the  EIC  star  detection  algorithms  in  an  effort  to  minimize  data 
loss.  Some  stars  detected  during  the  very  early  stages  of  EIC  development 
are  nc  longer  detected  yet  they  have  been  confirmed  as  real  by  ground-based 
telescope  observation. 
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COLUMN  HEADINGS  FOR  EIC-1 


EIC 


RA, 

DEC 


FL 


SDRA, 

SDDC, 

SDFL 


OBS 


1234 


Star  Number 

EIC-1  star  numbers  are  assigned  in  right  ascension  order  for 
1950.  0  coordinates. 

Coordinates 

Measured  right  ascensions  and  declinations  of  the  EIC-1  sources, 
corrected  as  described  in  the  text  and  precessed  to  equinox  1950.0. 
Right  ascension  is  listed  in  hours,  minutes,  seconds  of  time  and 
declination  is  listed  in  degrees,  minutes,  seconds  of  arc. 

Flux  Density 

The  flux  density  is  the  average  of  the  2.  7-pm  representative  mea¬ 
surements  in  units  of  10  *  W  cm  ^  pm  *. 

Standard  Deviations 

The  sample  standard  deviations  are  listed  for  RA,  DC,  and  FL. 

The  right  ascension  standard  deviation  is  listed  in  seconds  of  time; 

the  declination  standard  deviation,  in  seconds  of  arc;  and  the  flux 

- 1  _2  - 1 

density  standard  deviation,  in  units  of  10  W  cm  pm 
Number  of  Observations 

The  number  of  observations  used  in  computing  the  tabulated  values 
of  right  ascension,  declination,  and  flux  density. 

Observation  Interval  and  Flux 

Each  column  corresponds  to  an  observation  interval  given  below. 

A  means  no  observations  for  that  interval.  Other  symbols 
represent  the  infrared  flux  in  10  percent  steps  of  the  average  value 
(given under  the  heading  FL).  For  example,  "A,"  the  tenth  symbol, 
means  (1.  0  FL)  ±  5%,  "5"  means  (0.  5  FL)  ±  5%,  "P,"  the  25th 
symbol,  means  (2.  5  FL)  ±  5%,  and  so  forth.  For  reference, 
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The  means  greater  than  (3.  5  FL)  ±  5%. 


Obse  rvation 
Interval 

Start  Date 

Stop  Date 

Observation 
Time  (hr) 

1 

9  Nov  1976 

10  Nov  1976 

14.  2 

2 

18  Feb  1977 

19  Feb  1977 

40.  5 

3 

21  Jul  1977 

23  Jul  1977 

69.9 

4 

16  Feb  1978 

16  Feb  1978 

23.  4 

EAG, 

ESA, 

ETM, 

EAF 


TYPE 


L,  B 


The  distance  between  EIC  coordinates  and  AGK3  coordinates  is 
listed  under  the  heading  EAG.  Similarly,  SAO,  TMSS,  and  AFGL 
coordinate  differences  are  listed  under  the  headings  ESA,  EIR, 
and  EAF,  respectively.  The  distance  is  the  great  circle  arc  dis¬ 
tance,  given  in  arc  sec,  EIC  minus  the  other  position. 

The  value  for  spectral  type  is  obtained  as  described  in  the  text  from 
HD,  AGIO,  or  SAO. 

Galactic  Coordinates 

The  galactic  longitude  (L)  and  the  galactic  latitude  (B)  in  degrees 
are  computed  based  on  the  following  definitions  for  right  ascension 
and  declination  of  the  galactic  north  pole  (chq,  6q)  and  the  galactic 
longitude  of  the  celestial  north  pole  (Lq)  for  a  given  year  (YEAR). 

aG  =  12h49m  +  8^  13  x  10‘4  (YEAR-1950) 

6 G  =  27!  4  -  5!  44  x  10‘3  (YEAR-1950) 

Lg  =  123°  -l!  33  x  10-3  (YEAR-1950) 


SAO,  Associated  star  numbers  from  the  SAO,  AGK3,  TMSS,  AFGL,  or  HD. 

AGK3, 

TMSS, 

AFGL, 

HD 

DM  Durchmusterung  indication,  zone,  and  number.  Indication  is 

abbreviated  as  BD  for  Bonner,  CD  for  Cordoba,  and  CP  for  Cape 
Photographic.  DM  data  are  obtained  from  the  HD-SAO-DM  Cross 
Index  or  the  AGIO,  in  that  order  of  preference. 
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